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Two >130-meter-diameter impact craters formed on Mars during the later half of 2021. These are the 
two largest fresh impact craters discovered by the Mars Reconnaissance Orbiter since operations 
started 16 years ago. The impacts created two of the largest seismic events (magnitudes greater than 4) 
recorded by InSight during its 3-year mission. The combination of orbital imagery and seismic ground 
motion enables the investigation of subsurface and atmospheric energy partitioning of the impact 
process on a planet with a thin atmosphere and the first direct test of martian deep-interior seismic 
models with known event distances. The impact at 35°N excavated blocks of water ice, which is the 


lowest latitude at which ice has been directly observed on Mars. 


eismic recordings of hypervelocity im- 

pacts (>3 km/s) are rare despite being 

the most common terrain modification 

process in the Solar System. Earth is 

shielded by its atmosphere, consequent- 
ly there are few seismically recorded ground 
impacts, and meteoroids that do reach the 
ground usually travel at terminal subsonic 
velocity and only form small craters (7-3). The 
Apollo Passive Seismic Experiments on the 
Moon recorded ground motions from artifi- 
cial impacts, but these had slow relative velo- 
cities (<2.6 km/s), with respect to typical 
impact velocities of comets or asteroids collid- 
ing with the Moon, and formed craters smaller 
than 30 m in diameter (4). Larger natural im- 
pacts on the Moon were detected but have not 
been associated with imaged craters (4, 5), and 
all are expected to be smaller than 100 m in 
diameter (6). On Earth, a multitude of seis- 
mic events with known source locations, for 
example, explosion sources, have been used 
extensively for evaluating seismic velocity 
models, even down to the Earth’s core (7). In 
contrast, there have been only a few confirmed 


seismic source locations on Mars (all impacts), 
but these were small (<12 m in diameter) and 
near InSight (<300 km away), so the seismic 
paths only sampled the shallow crust (8). The 
two newly formed impact craters reported here 
allow for an evaluation of deep interior Mars 
global velocity models and observations of the 
dynamics of the hypervelocity impact process. 

The notable impacts (Fig. 1) were discovered 
using the Mars Reconnaissance Orbiter (MRO) 
Context Camera (CTX) (9) and the Seismic Ex- 
periment for Interior Structure (SEIS) (0) of 
the Interior Exploration using Seismic Investi- 
gations, Geodesy and Heat Transport (InSight) 
mission (17). Our CTX team independently dis- 
covered the Amazonis crater that we associated 
with the S1094b event. The earlier S1000a 
event, given its similar seismic signature, was 
then used to direct a search with MRO cam- 
eras to find the Tempe impact crater site. Both 
impacts generated craters >130 m in diameter, 
making them the largest fresh craters identi- 
fied since the beginning of the MRO mission 
16 years ago. The seismic events have identi- 
fiable surface waves, distinguishing them from 


other recorded and analyzed events on Mars 
and indicating shallow sources (12). Before 
these events, surface waves had not been un- 
ambiguously identified on any terrestrial planet 
other than Earth. The closer impact (S1094b) 
occurred at a distance of 58.5° (3460 km) from 
the InSight lander on 24 December 2021 and 
formed the larger of the two craters (150 + 10 m 
in diameter). The other impact (S1000a) oc- 
curred at a distance of 126° (7455 km) from the 
InSight lander on 18 September 2021 and 
formed a cluster of craters (the largest being 
130 + 12 m in diameter). The formation of the 
craters was time-constrained using the MRO 
Mars Color Imager (MARCI) (13) to within a 
day (table S1), making the association with the 
seismic events highly probable. The seismic 
events associated with the impacts have sim- 
ilar characteristics, both with ~4.0 magnitudes 
(tables S1 and S2). Because the seismic waves 
traveled deep in the mantle, both events are 
critical for analysis of mantle velocity models. 
However, given the S1000a event’s long dis- 
tance from the lander, direct seismic waves are 
eclipsed by Mars’ core (14) and more-complex 
bouncing seismic body wave phases (PP and SS) 
were detected (fig. S1). The additional attenua- 
tion and scattering experienced by these waves 
obscure the source characteristics, making 
source analysis much more challenging. In 
addition, the S1000a-associated crater is lo- 
cated on the side of a graben (fig. S2), which 
perturbed the blast pattern and prevented an 
easy identification of impactor parameters. In 
contrast, the closer impact (S1094b) occurred 
in a flat, dust-covered region. We first analyze 
the impact process for this closer impact be- 
fore considering the implications for Mars in- 
terior models of both impacts. 

Prominent surface albedo disturbances 
surrounding the S1094b impact allow for the 
estimation of ephemeral events that would 
otherwise be unknown, such as the impactor 
trajectory and the extent of atmospheric blast 
waves (Fig. 2 and fig. S3). The bearing of the 
bolide was estimated to be 60° + 5° clockwise 
from north by measuring the up-range “for- 
bidden zone” (15) in the albedo ray pattern and 
a down-range extended cluster of secondary im- 
pacts. We infer that the impactor approached 
the surface at an elevation angle of ~30° from 
horizontal. A steeper angle requires the asym- 
metric ejecta pattern to be muted, and a much 
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Fig. 2. Orbital images of the impact crater and surrounding area. CTX 
image (main panel; image ID U05_073077_2154_XI_35N170W): The hyper- 
velocity impactor traveled from southwest to northeast at an inferred azimuth 
of ~60° (fig. S3), creating a Mach cone shock wave that altered the surface 


shallower angle should have led to an elliptical | one generated by the passage of the meteoroid 
crater planform (15, 16). Two distinct arcuate | through the atmosphere (Mach cone) and the 
rays (“scimitars”) extending approximately | other by the ground impact (17), thus indicat- 
northwest and south of the impact site likely | ing that both blast waves extended to at least 
formed from the superposition of two atmo- | 18 km laterally. These arcuate rays provide an 


spheric blast waves disturbing surface 
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dust: | independent, albeit consistent, measure of the 
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albedo up-range of the impact, region labeled A in the figure. The inner 
dark ring, near B, is interpreted to be the result of blast wave mobilization 
of surface fines, impact-derived material directly deposited on the surface, 
or by ejecta-induced disturbances of the surface dust. The absence of 


Fig. 1. Impact event location map. The location of 

the impact craters (diamonds) and the InSight lander 
(yellow triangle) are shown. The S1094b crater is 
located at 34.80°N, 189.92°E in Amazonis Planitia. 

The S1000a crater is located at 38.11°N, 280.12°E in 
Tempe Terra. The great-circle paths between the new 
craters and InSight are superimposed onto the underlying 
globe image derived from MARCI (13), Mars Orbiter 
Camera (49), and Mars Orbiter Laser Altimeter (MOLA) 
data (50). The seismic epicentral distance estimates 

are indicated by the dashed white lines that extend 

over the azimuthal uncertainty estimate. The inset 
shows MARCI images from before and after the impacts. 
The MARCI images have ~2 km per pixel resolution 

at nadir. 


up-range ejecta disturbances indicates 
an oblique (~30° elevation) impact (15). 
Faint arcuate rays, labeled C, emanating 
cross-track of the impactor were likely 
caused by the superposition of the 

Mach cone and the atmospheric blast 
(17), indicating that both blast waves 
propagated out at least 18 km. The long- 
range ejecta-induced disturbances are 
concentrated in the down-range direction, 
region D, extending to at least 37 km. HiRISE 
image (inset; image ID ESP_073077_2155): 
The crater has a rim-to-rim diameter 

of ~150 m. The crater floor has an 
irregular shape, with a depth of roughly 
21 m. The light-toned material, for example, 
areas indicated by arrows labeled E, 
around the crater is inferred to be water 
ice ejected during the impact. 
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impactor trajectory (56°). The meteoroid struck 
the surface at 18:49 LMST (Local mean solar 
time, thus impacting on the orbital trailing 
side of Mars. We estimate the radial extent of 
surface dust disturbance from the crater to be 
9 km (fig. S3). This limit is consistent with the 
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Fig. 3. Seismic observation of S1094b, using de-glitched broadband data. (A) Vertical component 
velocity spectrogram. The event occurred at the end of a noisy period typical of martian afternoons. 

(B) Vertical component velocity time series bandpassed between 1 and 10 s and the derived spectral 
envelope. Phase picks for P, S, and Rayleigh (R1) wave arrivals are indicated with pick uncertainty indicated 
by black bars in the time series. (C) Waveform details of the P and S body waves (left) and Rayleigh 

wave (right). (D) Displacement spectra for the P, S, and Rayleigh waves and pre-event noise. See fig. S1 for 


a similar analysis of S1000a. 


atmospheric blast pressure produced by a 0.1 
to 1 kiloton (4 x 10™ to 4 x 10” J) surface 
explosive source (fig. S4). The size of the atmo- 
spheric blast allows us to evaluate its contri- 
bution to the seismic signal. 

Crater size is an important quantity for 
estimating the kinetic energy and momentum 
of the impactor for use in numerical models. 
An image from another camera on MRO, the 
High Resolution Imaging Science Experiment 
(HiRISE) (18), revealed additional details of 
the crater and its immediate surroundings 
(Fig. 2). The crater is irregular in shape, with 
an estimated rim-to-rim diameter of 150 + 10 m. 
Its depth, measured from crater floor to crater 
rim on the basis of photogrammetry results 
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using HiRISE stereo images, is roughly 21 m. 
The abundant craters surrounding this impact 
are likely almost all secondaries generated by 
the primary impactor, as, in comparison, areas 
far (>10 to 20 km) from the new impact have 
few small craters. The bright patches and blocks 
surrounding the crater reveal that it excavated 
water ice from the subsurface at a lower lati- 
tude (35°N) than any prior ice-exposing crater 
(39°N) (19). 

The geological context from orbital imagery 
aids in determining the appropriate physical 
models to use in numerical calculations. The 
$1094b crater is in the Amazonis Planitia re- 
gion in an area of rugged volcanic plains (20), 
with CTX images showing lava-flow morphol- 
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ogies mantled by a modest cover of debris. 
The lava flow indicates that a target ground 
with properties of porous fractured basalt is 
appropriate for modeling the surface impact. 
To account for a harder rock site at the crater 
compared with the region around InSight, we 
use a local subsurface velocity model based on 
terrestrial lava flows (27) extrapolated to Mars 
surface conditions (22, 23) (figs. S5 and S6) for 
seismic modeling. 

The seismic source duration for an impact 
of this size is expected to be shorter than the 
crater formation time scale and limited to the 
duration of the nonlinear shock wave propaga- 
tion regime (6). The seismic event S1094b had 
a very broad frequency content and relatively 
flat spectra, from 0.1 Hz to 3 Hz (Fig. 3), with a 
signal lasting >100 min (24) owing to propaga- 
tion coda (25). The event had an impulsive 
first-arrival P wave (1-s uncertainty), followed 
by an emergent strong S wave 6 min later 
(20-s uncertainty). The third wave observed, 
arriving 8 min later, was a Rayleigh surface 
wave, expressed as a long-period dispersed 
pulse with an 8- to 15-s period (72). All body 
wave phases are characterized by a long coda, 
indicative of strong scattering due to a near- 
surface source. The spectra display unusually 
high corner frequencies (3 Hz) compared with 
most other seismic events recorded on Mars. 
Shock physics modeling of the impact in a 
porous fractured basalt target (Fig. 4) indicates 
that most of the seismic moment was con- 
tained within a few hundred meters of the im- 
pact (Fig. 4A, orange bar). The moment release 
occurred over a short time period consistent 
with the cutoff frequency of 3 Hz identified in 
the P-wave displacement spectrum (Fig. 3D). 
Above the cutoff frequency, the amplitude 
shows a cubed frequency drop-off, as also ob- 
served for closer, smaller impacts on Mars (8) 
and for shallow explosions on Earth (26, 27). 

Seismic moment (Mo) is the key quantity 
that links the orbital observations of the im- 
pact and the impactor parameters to the seis- 
mic observations. For $1094b, the seismic 
moment estimate from S body waves was 
originally calculated assuming a marsquake at 
50 + 30 km depth to be 1.3 x 10” N-m (28). 
However, the impact seismic source deposits 
its energy at much shallower depths, in the 
strongly shocked region, estimated by impact 
modeling to be at a depth of between 17 and 
120 m, or ~50 m (Fig. 4B). We use our lava-flow 
seismic velocity model to conclude that the 
moment of a source at this depth is ~100 times 
smaller than a corresponding deep-crustal 
source for the same observed amplitude (Fig. 
4A and fig. S5). This is comparable to the seis- 
mic moment estimated from surface wave spec- 
tra, Mo = 7.5 x 10° N-m, at the same source 
depth (fig. S7). 

Empirical and numerical models were used 
to compute seismic moments for S1094b on 
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meteoroid into fractured basalt. Two snapshots, at 50 a 


nd 160 ms, show the zone of seismic wave generation 


where the shock pressure (Pshock) is substantially higher than the lithostatic pressure (Pith). 


the basis of the observed crater size. The im- 
aged crater diameter of 150 + 10 m corre- 
sponds to a vertical impactor momentum of 
3.3 + 1.4 x 10° N-s according to empirical 
crater-scaling relationships and impactor 
mass, angle, and velocity probability distribu- 
tions (figs. S8 and S9). For these values, nu- 
merical simulations predict seismic moments 
of 0.5 x 10"? to 1.2 x 10’? N-m for impacts in 
regolith and 2.8 x 10° to 7 x 10° N-m for 
fractured rock conditions (29, 30). These 
estimates are consistent with the observed 
seismic moment corrected for relevant depths 
in our subsurface model (Fig. 4A). The seis- 
mic efficiency was estimated to be 10 on 
the basis of scaling relations between seismic 
moment and crater diameter (37) with an order 
of magnitude uncertainty. This is lower than 
values estimated for lunar and Earth analogs 
(32) but larger than that previously modeled 
for small martian impact craters (29). 
The extensive blast pattern around the 
$1094b crater suggests that some of the seis- 
mic energy may have also originated from en- 
ergy released in the atmosphere and then 
coupled to the ground. Numerical impact sim- 
ulations suggest that up to 10% of the impact 
energy was partitioned into kinetic energy on 


(33). For impact scenarios on Mars similar to 
$1094b, simulations with an atmosphere fur- 
ther suggest that ~5% of the impact energy is 
partitioned into the blast wave (34). An esti- 
mated impact energy of between 1 x 10” and 
8 x 10” J (fig. S4) could produce an atmo- 
spheric blast comparable to a 0.1 to 1 kiloton 
(4 x 10" to 4 x 10° J) surface explosion. 
Therefore, both seismic and image observa- 
tions are consistent with such a blast and pro- 
vide coherent constraints in time and space, 
respectively. Semiempirical airblast theory (35) 
extrapolated to Mars suggests that such a blast 
would transition from the strong shock regime 
after 0.2 to 0.4 s, which is consistent with the 
observed ~3 Hz P-wave cut-off frequency. The 
induced blast pressure is sufficient to mobilize 
surface dust to a radius of ~10 km, which is 
consistent with the observed disturbed dust 
pattern (Fig. 1 and figs. S3 and S4). The esti- 
mated blast energy translates to an atmo- 
spheric moment Mg = (y — DE = 0.1 x 10” to 
1.3 x 10” N-m, where y is the adiabatic index 
(1.33 for Mars), and F is the blast energy (36). 
This is remarkably consistent with the moment- 
depth model when extrapolated to the sur- 
face (Fig. 4A), implying that a moment of only 
10° N-m released in the atmosphere could 


the planetary surface, primarily in the ejecta 
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explain a substantial part of the seismic body 
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wave signal, with the remaining part coming 
from direct coupling of the impactor with the 
ground. 

On Earth, atmospheric explosions easily ex- 
cite surface waves (37) and are highly sensitive 
to burst altitude (38). Solid Mars Rayleigh 
modes with atmospheric coupling (74, 39) are 
predicted to have excitation coefficients up to 
10 times larger for a near-surface atmospheric 
source compared with one that is 50 m below 
the subsurface (fig. S7). Surface waves are ex- 
pected to have an increase in the excitation 
coefficients between 0.1 and 0.15 Hz for sources 
above 50 m altitude. The increase is not ob- 
served in the estimated S1094b spectra and 
may be due to attenuation and scattering, 
which is not unexpected given that scattering 
effects are predicted to generate increasing 
attenuation of surface waves with frequency 
on the Moon and Mars (40). Comparison of 
$1094b surface wave spectra with near-surface 
excitation coefficients (fig. S7) suggests that 
a portion of the surface wave signal could 
have originated from the blast just above the 
surface. 

The Marsquake Service (MQS) (41) using 
SEIS data (24) estimated the seismic locations 
for both events (S1094b and S1000a). The dis- 
tance to the events is determined using 
S-minus-P arrival times (SS-minus-PP for 
$1000a) (42, 43) and polarization measure- 
ments of P and Rayleigh waves (PP for $1000a) 
(72). For the closer event, S1094b, the epi- 
central distance from the InSight lander was 
estimated to be 59.7° + 6.1° (3530 + 360 km), 
as compared with the actual distance of 58.5° 
(3460 km), a difference of only 70 km. For the 
second impact $1000a, the distance was esti- 
mated at 128.3° + 19° (7591 + 1240 km), as 
compared with the actual distance of 126.19, 
a difference of 130 km. Additional source 
parameters for these events are detailed in 
tables S1 and S2. 

The close agreement between distance esti- 
mates and the imaged locations increases our 
confidence in the martian seismic velocity 
models (44-48) for the regions sampled by 
the direct body waves (fig. S10). In particular, 
the models indicate the absence of mantle 
discontinuities in the 600 to 700 km depth 
range, which is the depth at which the P, S, 
PP, and SS waves turn (44). For the $1000a 
event, the Paige phase, the P wave that diffracts 
along the core mantle boundary (CMB), has 
been tentatively identified (74). The PP-Paige 
travel time difference is sensitive to lower 
mantle P velocities below 800 km. Current 
models at these depths are constrained by core- 
reflected S phases and the mineral physics- 
based Vp/Vs ratio (fig. S10). The observed 
PP-Pai¢g does not match the predicted values 
given by these models. This mismatch implies 
that either the P velocities at the CMB need ad- 
justment, or the Vp/Vs ratio in the lower mantle 
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is different from current predictions. These two 
events act as calibrated measurements and 
help select among various martian interior 
seismic velocity models (44-48); they corrob- 
orate Mars mantle velocity models to 800 km 
depth and will help to improve future models 
down to the CMB. 


The first two recorded teleseismic events 


on Mars with orbital ground-truth observa- 
tions have been used to constrain martian 
interior seismic velocity models and infer 
dynamic impact processes including seismic 
moment release, impact source duration, and 
atmosphere-subsurface energy partitioning. 
The success in observing the formation of 
impact craters on Mars using instruments 
on several missions opens up a more detailed 
understanding of impact dynamics, atmo- 
spheric physics, and the exploration of plan- 
etary interiors. 
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An insightful impact 

On 24 December 2021, the seismometer for the InSight mission on Mars detected a large seismic event with a distinct 
signature. Posiolova et al. discovered that the event was caused by a meteor impact on the surface of Mars, which was 
confirmed by satellite observations of a newly formed 150-kilometer crater. The surface nature and size of the impact 
allowed Kim et al. to detect surface waves from the event, which have yet to be observed on Mars. These surface 
waves help to untangle the structure of the Martian crust, which has various amounts of volcanic and sedimentary 
rock, along with subsurface ice, in different regions of the planet (see the Perspective by Yang and Chen). The 
characteristics of the impact itself are important because they provide a seismic fingerprint of an impact event that is 
different from the marsquakes observed so far. —BG 
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